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Young's modulus, as a function of temperature up to 1000 °C, was determined for eight 
rutile specimens of different crystallographic orientations. From these, the following four 
elastic compliances or combinations of them were determined as a function of temperature, 
Su, Sxi, 2.s'i:{-l-.s'44, and -sn — ,si2 — .s'fi(i/2. Tlui method used was the same as had been used in 
a previous study for rutile at room temperature. 



1. Introduction 

In a previous investigation [1],^ all the elastic con- 
stants of rutile were determined at room temperature 
by an application of the resonance technique. 
Briefly, the method of approach consists in deter- 
mining the mechanical resonance frequencies of both 
flexural and torsional vibration of a group of cylin- 
drical rods of different orientations of the same type 
of crystal. Young's and the shear modulus are com- 
puted for each bar from these frequencies and from 
its dimensions and density. Then, from a set of 
equations given by Nye [2], for histance, relating 
Young's and the shear modulus to orientation and 
elastic compliances, one solves for the elastic com- 
pliances (since the other parameters. Young's modu- 
lus, shear modulus and orientation are known) . 

The number of equations required in the set, and, 
hence, also the number of specimens of different 
orientation, depends upon the number of elastic 
constants required to completely describe a particular 
crystal. If more than the minimum number of 
specimens are available, then, from an overdeter- 
mined set of equations, one may obtain an estimate 
of the precision of measurement of the elastic con- 
stants as well as the elastic constants themselves. 
This has been done for corundum [3] as well as 
rutile [1]. 

In this paper, some of the elastic compliances 
already determined at room temperature for rutile — 
those arising from the flexural resonance frequen- 
cies — are extended to elevated temperatures using 
the same method of approach. As far as we know, 
there is no data available on the temperature depend- 
ence of any of the elastic properties of rutile. 



1 Figures in brackets indicate the literature references at the end of this 
paper. 



2. Experimental Procedure 

The method for obtaining resonance frequencies, 
both at room and elevated temperatures, has already 
been described [4-6]. However, certain refinements 
(such as the use of fine wires of phosphor bronze as 
suspexision members) which were developed to obtain 
torsional resonance for the shorter rutile specimens 
[1] were not applicable at elevated temperatures. 
Furthermore, with the glass and wire fibers used at 
elevated temperatures as suspension members in 
place of the cotton or silk threads used at room 
temperature, it was not possible (maiidy because of 
the small cross section of the specimens) to obtain 
satisfactory torsional resonances even for the longer 
specimens. Consequently, this investigation is re- 
stricted to the flexural resonance frequencies and to 
those elastic compliances which can be derived from 
them, as mentioned above. 

For this purpose, 8 of the original 16 specimens 
were used. This included aU of the longer ones 
(which yielded the most reliable values of resonance 
frequency) and those of the shorter ones having 
orientations which would be most significant in giving 
a Young's modulus-orientation relationship. The 
specimens used, along with their length, nniss, and 
orientation, are given in table 1. The same notation 
for designating specimens and orientation, as well as 
most of the other notation used in the previous 
paper [1] is retained here. 

Flexural resonance frequencies for each specimen 
were determined in steps of about 100 °C, allowing 
sufficient time for thermal equilibrium to be attained 
at each temperature. Usually two or three sets of 
temperature determinations were taken for each 
specimen. The glass fiber suspensions (which yielded 
the more reliable results) were used for all the speci- 
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mens up to the limit of their appHcabihty; from 700 
to 800 °C. The wire suspensions gave satisfactory 
responses with the longer specimens up to the highest 
temperatures of the furnace, about 1400 °C, but were 
unsatisfactorv for two of the sm.aller specimens over 
about 800 °C. Since the changes in elastic compli- 
ances were fairly gradual and continuous with 
temperature, it was felt to be a more accurate pro- 
cedure to extrapolate the values obtained with the 
glass suspension for these two specimens up to 1000 
°C rather than use the data obtained with the wire 
suspension. 



Table 1. Properties of ^^ single crysiaV^ rutile speci?nens 





Mass 


Length 


End 


Center 


other end 


Speci- 
men 
number 


d 


</> 


d 


6 


4> 




All angles in degrees 



Linde rods 



23 

24 



(1 

2. 6317 


cm 
13. 236 


/ 12.5 
I 12.3 


2.6 
2.2 


|l2.7 


2.1 


f 14.3 
1 14.2 


2. 8205 


12. 168 


14.0 


7.3 


16.1 


7.0 


14.8 


2. 3654 


11.714 


43.3 


34.1 


45.6 


37.9 


44.7 


1. 1523 


5.816 


58.7 


24.0 


58.5 


24.2 


59.3 


15. 0251 


9.848 


88.8 


44.1 


88.4 


43.4 


90.0 



1.8 
1.5 
6.3 
33.6 
23.7 
44.5 



NBS rods 


49 
50 
51 


7. 1208 
21. 5907 
16. 0263 


11. 128 
17. 107 
17. 634 


87.8 
88.8 
86.2 


0.0 
5.4 
0.4 


87.8 
89.5 
86.2 


0.0 
5.4 
1.0 


88.0 
90.0 

86.8 


1.5 
5.0 
0.9 



Young's modules at some elevated temperature, 
Ytj was computed from the resonance frequencies, 
from the equation 



Yt= Yo{ft/Joy{lt/lor{mo/m^) (no/n,y 



(1) 



The room temperature value of Young's modulus, 
Fo, is known from the previous investigation; 
ft and /o are the flexural resonance frequencies at tem- 
perature t and room temperature respectively. It 
should be noted that the ratio, Jt/fo, applies to over- 
tones as well as to the fundamental resonance 
frequency. This was used for the longer specimens 
to provide additional points (obtained from over- 
tones as well as the fundamental) during a modulus- 
temperature run. The points based on overtones 
also served as an internal check on the consistency 
of the modulus-temperature relationship. The 
terms involving I, m, and n are the correction for 
thermal expansion. Since the thermal expansion 
of rutile is anisotropic, the coefficient being about 
30 percent larger in the ^^c" or [001] direction than 
in the '^a" or [100] direction, it was necessary to take 
this into account in a manner which has already 
been described [7]. The length of the long axis of the 
specimen at temperature t is It, and Iq is the same 
dimension at room temperature; similarly for m 
and n, n being the cross-sectional dimension in the 



direction of flexural vibration and m being the cross- 
sectional dimension perpendicular to n. A more 
detailed description of the application of this aniso- 
tropic thermal expansion correction to rutile is given 
in the appendix. Actual values of the thermal 
expansion of rutile in both '^(i" and ^^c" directions 
were supplied by Richard K. Kirby of the NBS 
staff. The possibility that certain orientations 
might change with increasing temperature due to 
this anisotropic thermal expansion was also con- 
sidered. However, it was found that even in the 
extreme case (^=1000 °C and ^=45 °C) the change 
in orientation would be insignificant, i.e., less 
than K^ 

The scatter of the experimental points of Yt as a 
function of temperature, based on repeated runs 
and (where used) on overtones, about a smooth 
curve drawn through the points was estimated to be 
about 5 parts per 1000. 



3. Results 

Values of 1/F^ interpolated from the experimental 
data at even 100 °C temperature intervals for the 8 
specimens measured are shown in table 2. The data 
at 25 °C are from the previous investigation [1]. 
The compliances are obtained from eq (6) of [1], 

^lYt=s^z=Su sin^ ^+533 cos^ d-\- (2^13+544) sin^ 6 cos^ d 

+1(^12-^11+^) sin^^sin^ 20. (2) 

This equation was solved by computer at each 
temperature (100 °C, 200 "^C. . . 1000 °C) in the 
same manner as before. It is noted that since 8 
specimens are used and 4 coefficients are solved for, 
an overdetermined set of equations is still available 
from which, as previously mentioned, one may 
derive an estimate of the precision of the coefficients 
from a least square best fit, as well as the coefficients 
themselves. Table 3 presents the four elastic com- 
pliances along with their standard deviations, com- 
puted in this manner. The standard deviations of 
the compliances, especially at the lower tempera- 
tures, compare favorably with those computed at 
room temperature. As might be expected, these 
standard deviations tend to grow somewhat larger 
at the higher temperatures. The four coefficients 
are shown graphically as a function of temperature 
in figure 1. 

Once the compliances are known, they may be 
resubstituted in equation (6) and ^theoretical" or 
^^computed" values of 1/F^ may be solved for at 
each temperature. The numbers associated with 
each value of 1/F^ in table 2 represents the difference 
between these ^ ^computed" values and the ' ^experi- 
mental" or '^observed" values given in the table. 
Following the same convention as before, a -\- sign 
indicates that the experimental value of 1/F^ is 
larger than the theoretical one and vice versa. 
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Table 2. Reciprocal of Young's modulus for rutile specimens of various orientations as a function of temperature, 

HYt=sl^ in 10-i2m2/iV (10-i3cm2/dyne) at indicated temperature (°C) 



Specimen No. 


25 


100 


200 


300 


400 


500 


600 


700 


800 


900 


1000 


21 


2.666 
2.690 
3.047 
3.881 
2.694 
6.777 
6. 650 
6.757 


2.708 
. 0005* 

2. 731 
-.0006 

3.092 
.0003 

3. 921 

-.0003 

2.778 
.0001 

6.777 
.0016 

6. 651 
-.0010 

6.755 
-.0006 


2.762 
.0015 

2. 783 
-.0012 

3. 147 
-.0013 

3. 973 
.0017 

2.872 
-.0002 

6.779 
. 0009 

6.656 
-.0017 

6. 759 
.0004 


2.815 
-.0001 

2. 840 
-.0013 

3. 203 
- . 0056 

4.032 
-.0070 

2.960 
-.0006 

6. 798 
. 0093 

6. 665 
-.0057 

6.764 
-.0054 


2.871 
-.0003 

2. 895 
. 0009 

3. 257 
-.0029 

4.074 
.0035 

3.048 
-.0004 

6.800 
.0022 

6. 681 
-.0014 

6.777 
-.0017 


2.927 
-.0001 

2.950 
.0005 

3.314 
-.0020 

4. 125 
.0024 

3.135 
-.0003 

6.820 
.0017 

6.704 
-.0010 

6.798 
-.0013 


2.986 
-.0016 

3.011 
.0016 

3.373 
.0004 

4.175 
-.0006 

3.223 
.0001 

6.846 
.0016 

6.732 
-.0010 

6.825 
-.0005 


3.050 
-.0011 

3. 072 
-.0005 

3.444 
.0081 

4.226 
-.0098 

3.320 
.0011 

6.877 
.0014 

6.769 
.0027 

6. 855 
-.0018 


3.108 
-.0023 

3.134 
.0019 

3. 512 
. 0026 

4.303 
-.0032 

3.408 
.0003 

6.922 
.0058 

6.812 
.0036 

6.889 
-.0086 


3.176 
-.0008 

3.198 
-.0002 

3.580 
.0054 

4.364 
-.0065 

3.501 
.0007 

6.964 
.0006 

6.859 
.0018 

6.944 
-. 0009 


3.241 


23 


-.0016 
3. 263 


24 


-.0004 
3. 647 


30 -- 


.0105 
4.423 


44 - 


-.0128 
3.597 


49 


.0019 
7.047 


50 


.0212 
6.911 


51 


-.0096 
6.998 




-.0092 



*Lower number in each box is difference between "observed" value of 1/Yt and "computed" value. Observed values of 1/Yt are given in all cases except at 25° C 
where computed values from pr(>vious investigation [1] are given. 



Table 3. Elastic compliances of rutile as a function of 
temperature in 10~^^ m^/N {10~^^ cm^jdyne) 



Temperature 


Sn 


S33 


2Sn-\-Sl^ 


Sn Si2 2 


°C 
25 


6. 787 ±. 0003* 
6.786 ±.0006 
6.789 ±.0010 
6.799 ±.0040 
6.808 ±.0016 
6.828 ±.0012 
6.855 ±.0009 
6. 886 ±. 0039 
6.926 ±.0035 
6. 973 ±. 0026 
7.036 ±.0088 


2.592 ±.0004 
2.635 ±.0010 
2.689 ±.0015 

2. 744 ±. 0066 
2.803 ±.0024 
2.860 ±.0017 
2.923 ±.0013 
2.988 ±.0058 

3. 045 ±. 0052 
3. 113 ±.0038 
3.181 ±.0131 


6.464 ±.0016 
6.528 ±.0039 
6.617 ±.0063 
6. 725 ±. 0273 
6.793 ±.0099 
6.881 ±.0071 

6. 965 ±. 005r) 

7. 065 ±. 0238 
7. 219 ±. 0217 
7.325 ±.0157 
7. 414 ±. 0540 


8.197 ±.0011 


100 


8.026 ±.0056 


200 

300 


7. 843 ±. 0041 
7.687 ±.0178 


400 


7.529 ±.0065 


500 


7.396 ±.0046 


600- 


7. 272 ±. 0036 


700 


7.142 ±.0155 


800 

900 


7.046 ±.0141 
6.955 ±.0103 


1000 


6. 890 ±. 0353 



•^ standard deviation. 




100 200 300 400 500 600 700 800 900 1000 

TEMPERATURE , °C 

Figure 1. Elastic compliance of rutile as a function of 
temperature. 



Figure 2 shows the Young's inodulus-orientation 
dependence for rutile at two selected temperatures, 
400 °C and 1000 °C, obtained in the same manner 



as the room temperature curves, which are repro- 
duced from the original study for comparison. 

It may be seen from the figure that in addition to 
the expected decrease in modulus for all orientations 
with increasing temperature, tlie anisotropy of the 
crystal decreases also. For any given value of d, 
the change in Young's modulus, as a function of (^, is 
less at 1000 ""C than at 25 °C. Also, for any par- 
ticular value of <^, the variation in 6 appears to 
decrease at the higher temperatures, i.e., the lines 
become flatter. 
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Figure 2. Youncfs modulus of rutile as a function of orientation 
at selected temperatures. 

e is angle of specimen axis with respect to "c" axis of crystal; 4> is angle between 
the projection of specimen axis in (001) plane of crystal, and "a" axis of crystal. 
(The subscript "/" is retained in Yf from the previous investigation [1] to indicate 
the "free" rather than the "pure" modulus.) 
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4. Appendix 

The problem is to express eq (1) in terms of the 
known thermal expansions for rutile for given 
crystallographic directions. We define these two 
thermal expansions for the ^^c'^ and ''a'' directions 
respectively as LtdLoc and LJLoa- 

(1) We consider first a ^'c'' axis rod, i.e., a speci- 
men of which the long dimension is in the direction 



of the crystallographic 



^/^^^and- 



mo 



For this condition, 

_nt__Lta 



SO that equation (1) becomes 

' V/o/ \lJ \lJ ' 



Y —Y 



(2) We next consider an ''a" axis rod. In this 
case, since the thermal expansion of the cross- 
sectional dimensions are not equal, a specimen 
having a circular cross section at room temperature 
will develop a noncircular cross section as the 
temperature increases. This will give rise to two 
separate flexural resonance frequencies; one in the 
direction of the short cross-sectional dimension, the 
^^a" direction, and the other in the direction of the 
long cross-sectional dimension, the ^'c" direction. 
These two types of vibration are considered sepa- 
rately. For vibrations along the ''a" direction, 

h-^h^, E±=^, and ^=^ 



Lq IjQa ^0 -LiQ, 



riQ 



(since ^^n^^ is always the dimension in the direction 
of vibration). Equation (1) then becomes 

^'■=>'«(f(z7:> 

For vibrations in the "c" direction, 






'"■ ™i= Jls and ^= 



L 

■Litia ^0 






and eq (1) becomes 



Y, 






It is clear that any specimen having its long 
axis in the (001) plane (^==90°, ip from to 45°) will 
have the same thermal expansion corrections as the 
two just given, since the same considerations holding 
for an ^^a^^ axis rod would apply. 



It also follows from these considerations that rods 
having their long axis in the (001) plane should 
develop two noticeably different curves of (/^//o)^ as 
a function of temperature. This was observed 
experimentally; however, when the appropriate 
thermal expansion correction was applied to the 
points forming each curve (giving YtlY^ as a function 
of temperature), the two curves were brought into 
coincidence within experimental error. Conversely, 
for ^^c" axis rods, it was observed that even if two 
resonance frequencies were present at room tempera- 
ture, due to the fact that the cross section was not 
perfectly circular (at room temperature), never- 
theless, the two curves of {ft/foY as a function of 
temperature, based on these two room temperature 
resonance frequencies, coincided within experimental 
error. This was further verification of the validity 
of the thermal expansion corrections. 

(3) The final case to be considered is that of a rod 
having its long axis at some angle, 6, with the '^c" 
axis. In this case, since the "c^^ axis does not 
coincide with any specimen axis (/, m, or n) as has 
been true in the previous examples, the thermal 
expansions along /, m, or n will generally not equal 
Ltc/Lnc or Lta/Loa, but will lie between these values. 
The following equation [7] then applies. 



V/o/\ L,^\L, Lj 



[3 cos' 6>-sin2 6>(l+2 cos';//)] 



} 



where yp is the angle between the direction of vibra- 
tion and the plane containing the ^^c'^ axis and the 
long axis of the rod. If this equation is written for 
each of the two directions of vibration, and the 
average of these two equations is taken, the following 
relationship results: 



V/o/ave. i^ 






+(f-f)p— "'»'}■ 



This final equation is also applicable to the 
specialized orientations previously discussed. It is 
useful for the case in which the long axis of the rod 
coincides with the [001] direction or lies in the (001) 
plane, but in which the direction of vibration does 
not coincide with a crystallographic axis. 
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